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Radiation  of  a  Hertz  Dipole  Moving 
around  a  Conducting  Sphere 

John  Lam 


ABSTRACT 

An  exact  expression  is  obtained  for  the  field  of  a  horizontal 
Hertz  dipole  moving  around  a  conducting  sphere  in  a  circular  orbit. 
This  expression  is  in  the  form  of  a  double  series  which,  in  the  case 
of  a  large  sphere,  is  evaluated  by  the  method  of  Watson  transformation. 
The  geometrical  optics  field  is  separated  out  and  its  properties 
examined.  It  is  found  that  the  incident  and  reflected  waves  are  of 


different  frequencies. 


2. 


1.  Introduction 

The  problem  of  wave  scattering  by  a  spherical  object  occurs  in 
many  branches  of  science  and  engineering.  In  the  majority  of  cases  the 
source  of  the  vave  is  fixed  relative  to  the  sphere.  A  rigorous  solution 
is  in  general  obtainable  by  the  method  of  separation  of  variables, 
giving  the  scattered  vave  as  an  infinite  series  of  eigenfunctions.  If 
the  radius  of  the  sphere  is  large  compared  to  a  wavelength  of  the  inci¬ 
dent  vave,  the  convergence  of  the  series  is  vexy  slow.  The  method  of 
Watson  trans  format  ion  is  then  applied  to  convert  this  series  into  a  more 
rapidly  convergent  series  of  "residue  waves".  This  method  of  solution 
can  be  carried  over,  with  only  minor  extensions,  to  problems  involving 
moving  sources.  In  this  work  we  calculate  the  radiation  field  of  an 
oscillating  electric  dipole  revolving  around  a  perfectly  conducting 
sphere.  This  problem  is  an  idealization  of  the  situation  of  a  transmit¬ 
ting  antenna  carried  by' an  artificial  satellite  in  an  orbit  around  a 
planet . 

2.  The  Hertz  Potential 

We  set  up  a  spherical  polar  coordinate  system  whose  origin  coin¬ 
cides  with  the  center  of  the  conducting  sphere,  the  radius  of  the  sphere 
being  b  .  The  oscillating  electric  dipole  revolves  in  a  circular  orbit 
of  radius  a  in  the  x-y  plane,  given  by  r  ■  a  >  b  ,  6  *  n/2  .  To 
simplify  the  problem  we  assume  it  to  be  polarized  in  the  z-direction. 
Without  further  loss  of  generality  we  can  write  down  its  electric  polari¬ 


zation  as  follows: 


3 


— itu  t 

P(r,t)  -  ^ - —  6 ( r-a )  6(0  -  £)  6(0  -  Jit)  e  (2.1) 

r  sin  0  2  -z 


where  p  is  its  electric  dipole  moment*  fl  its  frequency  of  revolution 
around  the  sphere,  and  toQ  its  proper  frequency  of  oscillation  increased 
by  the  dilatation  factor  (l  -  fl2a^/c^)’^^. 

To  find  the  field  of  the  moving  dipole  ve  first  calculate  the  Hertz 
potential  n(r*t)  which  satisfies  the  wave  equation 


l  a2 

K-K)  £(r,t) 
c 


P(r.t) 


(2.2) 


Because  of  our  particular  choice  of  the  polarization  this  is  really  a 
scalar  equation  with 

n(r*t )  *  n(r,t)  e  (2.3) 

—  —  —  “Z 

The  angular  delta- function  6(0  -  fit)  is  a  periodic  function  and  has  the 
Fourier  series  representation 


4(0  -  at)  -  fe  j 


(2.4) 


Putting  (2.4)  in  (2.2)  and  expanding 


as  im0  -  i(u  +  mfl)t 

n(r,t)  ■  J  nm(r,0)  e 

m»-» 


(2.5) 


we  obtain  the  equation 


(  '  Sr2  '  *  r2sin  0  3«  ■i,‘  *  "  '  rCn2«  +  )  "“<r’S) 


- - 2 — X  6(r-a)  6(0  -  f) 

2  ire  a  * 
o 


(2.6) 


k 


For  convenience  ve  have  defined  the  quantities 

w  ■  w  ♦  mfl  ,  k  ■  k  +  mK 

no  no 

ko  "  wo/c  *  K  »  n/c  (2.7) 

To  solve  equation  (2.6)  ve  divide  the  whole  space  into  regions  1  and 
2  according  to  r<a  or  r  >  a  ,  respectively.  In  these  regions  the 
solutions  of  (2.6)  are 

"i(r*0)  "  J  Atm  htl)(kna)  Jl(kmr)  F?(cos  0)»  r  <  a 
t-0 

(2.8) 

nB<r»°)  "  l  Btn  Vkma)  hil)(knr)  Fi(c08  0)»  r  >  a 
t»0 

is  chosen  to  he  finite  at  the  origin,  and  n2  satisfies  the  radiation 
m  m 

condition  at  infinity.  The  constants  A.  and  are  determined  by 

1  2 

the  boundaxy  conditions' on  the  sphere  r*a  .  First  n  and  n  are 
continuous  across  the  boundary,  and  ve  get 


(2.9) 


Second  ve  integrate  both  sides  of  (2.2)  over  the  volume  of  a  flat  pill¬ 
box  of  infinitesimal  height  bounded  by  the  surfaces  r  *  a  -  e  and 
r  *  a  +  e  .  The  second  term  on  the  left-hand  side  gives  negligible  con¬ 
tributions  since  ll(r,t)  is  continuous  across  r  ■  a  .  Using  the  Gauss 
divergence  theorem  on  the  first  term,  ve  get 


5. 


J  V2n(r,t)  dV 

-  ||  a2sin  «  a»  d0  n2(r.t)|r_a+e  -  57  n1(r,t)lr.s.E]  (2.10) 


The  contribution  to  the  surface  integral  from  the  side  of  the  pill-box 
is  negligible.  On  the  other  hand  the  right-hand  side  of  (2.2)  gives 


_  P(r,t) 


dV 


-I! 


a  sin  0  dO  d0 


-iw  t 

j  0  6(0 


eo  a 


f)  6(0  -  Qt) 


(2.11) 


Comparing  (2.10)  and  (2.11)  we  have  the  boundary  condition 


3r 


r«a+e 


_3 

3r 


r»a-c 


-iw  t 

e  °  6( 9  -  £)  4(0  -  St) 


(2.12) 


The  right-hand  side  can  be  expanded  into  a  series  of  spherical  harmonics 
by  means  of  the  orthogonality  relation 


r  i  i(m-m' )0 

;  I  ;'(cos  0)  ,  (cos  0)  e  dfl 


Uir  (t  +  |m|  )l  .  . 

24+1  (l-  m  )T  °IV  °mm' 


(2.13) 


The  result  is  as  follows: 


-iw  t 


6(0  -  p  6(0  -  fit) 


“  1  /  ■  I  |  \,  im0-  iu  t 

i  i  *po)  *;(*.•). 


2  .A  L  6  ***  U+ 

c  a  1*0  m*-£, 
o 


(2.1U) 


Putting  (2.5),  (2.8),  (2.9)  and  (2.lfc)  into  (2.12)  we  get 


i &r  <2t+1> 


i+  m  TT  r 


ikm  • 


m|  *  l 


|m|  >  l 


Here  use  has  been  made  of  the  identity 


(2.15) 


jt(z)  h[l)'(z) 


J j(z )  h*l}(z) 


(2.16) 


Hence  we  have 


n1(r,t) 


n2(r,t) 


o  1*0  m*-l  1 


W>\  lm0-  U  t 

X  pJ(cos  0)  e  m  (2.17) 

'W>  hi1>(kmr>' 


We  notice  that  the  frequency  spectrum  of  n(ir,t)  consists  of  discrete  lines 


at  at  ■  u  ♦  ml)  . 
m  o 


7. 

3.  The  Incident  Fit Id 

The  incident  electric  and  magnetic  fields  are  derived  from  the  Hertz 
potential  according  to  the  equations 


Einc  -  7  x(7  x  n)  ,  cBinc  «  -  |r  V  x  n 
—  —  “  c  it 


(3.1) 


In  what  follows  we  only  need  to  know  the  radial  components.  These  are 
given  by 

.2 


cB 


inc  11  3 


r  c  3t  30 


n 


,inc 


2  cos  0  3  sin  0  3 

r  3r  2  30 

r 


["  8  In  »  3g 

["  r  Sr  98 

(.  _1 _ l  1  -£>1  r 

'  .18  6  39  8  39  .ln28 


cos  0 


(3.2) 


Substituting  (2.17)  into  (3-2)  we  have 


cB 


inc 


i  ®  1 

—  £  £  A,  ink 

r  .L.  L  .  £m  m 
£■1  mF-l 


hj1  (k  a)  J,(k  r)| 
t  In  l  a 

jt<v>  hi1)(v> 


P^(cos  0)e 


im0  -  iu it 
m 


4"° 


2  cos  0 


A  k 
OO  0 


h^l}(k  a)  JI(k  r) 
o  m  o  m 


JJM)  h!l)’(k  r) 


'om  o 


i  *  1 


r  £-1  m»-£  2£+1  »Jt(kma) 


-iu  t 
o 


(£-l)(£+|m|) 


Ji-iV’l 


+  (£+2)(£-|m|+l)  (  4+1  “  I  .(cos  0)1  e 


im0  -  iu  t 
m 


(3.3) 


i+l  m 


8. 


In  deriving  the  expression  for  ve  have  used  the  identities 

(1  -  z*)  F»'  (z) .  p^_i(2) .  »,IhL;K)  ^+1(.) 

2  A^l 

z  pPfz)  .  1+1*1  pP  (z)  ,  t-lmkl  J1  ,  * 

Z  *V  '  2£+l  *£-lU;  2C+1  *£+lU'  * 


i(a)  j 

c 

1  Jc-i(z)\ 

£  +  1  1 

^c+i^ 

"  2C+1 

■  2C+1 

J[(.)  , 

1 

j  J£-l^\ 

.  1 

1  Jc+i^z 

h[1)(Z)l 

"  2C+1 

2  C+l 

ki<a 

inc 

After  some  rearrangement  ve  can  write  E  in  the  simpler  form 

r 


.  i  08  £  /■'I'm-’  \ 

*  . 


J£(kmr)  t  im0  -  ini  t 

m 


h£l)(kmr) 


where 


C.  »  k 

Cm 


m  [AC-l,m 


(l+l)(t- |m | )  I ht-l^km°^ 1 
24-1 


(1), 


+  A. 


C(£+|m|+l)  /  h£+l^kma)  I] 


21  +  3  »WV> 


£+l,m 

Substituting  (2.15)  into  (3.6)  ve  get 

Ctm  *  V»0  -  lkm  [<fl>(U|»D  ^(0) 

*  /h^(ka)»i 

+  C(£-|m|+l)  Pj  (0)  ] 

lWkma>'J  ’ 


‘Sv 

Jt-l(kma) 


(3.U) 


(3.5) 


(3.6) 


(3.7) 


9. 


10. 


,inc 


r-e- 

Uire  r  , 

o  i, 


l  l  <?fu{£ 

-1  BF-t  K 


u- 

in 

III 

u+ 

m 

I7i 

ik 


•  —  t  A IV 

fp?  (o)-S  » 


i!  I 


W> 

h<l)(kmr) 


im0  -  ioj it 

P^(cos  0)  e  (3.11) 


In  the  second  equation  appears  to  be  discontinuous  on  the  sphere 

r  *  a  .  This  discontinuity  is  only  apparent.  The  difference  in  the 

inc 

egressions  for  E  on  both  sides  of  r  *  a  is  an  infinite  series 

r 

which  can  be  shown  to  converge  to  zero  everywhere  except  at  the  source. 
U.  The  Debye  Potentials 

To  describe  the  scattering  of  electromagnetic  waves  by  a  sphere  it 
is  most  convenient  to  work  with  the  Debye  potentials  u  and  v  ,  from 
which  the  fields  are  derived  according  to  the  formulas 

7  1  3u  v 

E*-V(r-—  ♦£  »  Tv) 

cB  -  -  Vx(r  x  Vu  -  r  i  |J)  (U.l) 

—  —  -  C  ot 


E  and  B  given  by  (U.l)  axe  solutions  of  Maxwell's  equations  in  free 
space  provided  that  u  and  v  separately  satisfy  the  scalar  wave  equation 


<’2  0 

c  3t 


Writing  out  (U.l)  in  component  form  we  get 


2  2 

,  3  1  3  \ 

rV 

3r  c  3t 


(U.2) 


11. 


1  1_  32(ru)  +  1  32(rv) 

r  sin  0  c  3t  30  r  3r  30 

1_  1_  32(ru)  1  32(rv) 

r  c  3t  30  r  sin  0  3r  30 


32  1  a2 

cBr  “  (“T  "  h  "TJ  ru 

3r  c  3t 

cB  .  I  + _ l _ I  sites! 

0  r  3r  30  r  sin  0  c  3t  30 


cB, 


1  32(ru)  1  1_  32(rv) 


0  r  sin  0  3r  30  r  c  3t  30 


(U.3) 


Thus  u  generates  a  magnetic  wave  and  v  an  electric  wave.  It  is  also 
clear  that  they  are  determined  by  the  radial  components  of  the  magnetic  and 
electric  fields  respectively. 

Since  u  and  v  are  solutions  of  the  scalar  wave  equation,  they 
have  the  eigenfunction  expansions 


os  i 

u  *  l  Id, 

A*1  m*-£  i 


'h^(k  a)  J,(k  r) 
l  m  ''i  m 


Jl(kma)  h[1>(1Imr) 


f^(co.  8)  e1"0-  lV 


I  I  * 

t«l  m*-t, 


Im 


[k  a  h^(k  a)],J.(k  r) 
m  l  m  m 

[k  a  J  (k  a)]’h|^(k  r) 
m  m  t  m 


im0  -  iw  t 

P?(cos  0)e 


(k.k) 


To  get  and  cBr  we  operate  on  ru  and  rv  with  the  operator 


2  2 
a  1  a 


2  2  2 
3r  c  3t 


as  indicated  in  (U.3).  Because  of  the  following  equation,  satisfied  by 


12. 


the  spherical  Bessel  and  Hankel  functions 


(4.5) 


this  operation  merely  brings  in  a  factor  i(fc+l)/r  .  Thus  comparing 
(3.11)  and  (4.4)  we  immediately  obtain  the  expressions  for  the  Debye 
potentials  of  the  incident  field: 


inc  P  “  r  21+1  U-  m  )l  ..2 


hl1)(kma)  W> 
Vkma)  h[1)(k»r) 


•  l 


im0  -  J  a)  t 

F^(cos  9)  e 


_inc  P  r  r  2i+l  (£-  m  )\  „m'  m 

y  iTTr  l  l  TTI+T)  Twin Ti  p»  (0)  ~ 


'o  £*1  m*-£ 


'V h,  V11'  W> 

!v  W*1'  hi1)(v> 


im0-  iu  t 
P?(cos  0)  e  m 


(4.6) 


i.  The  Scattered  Field 


The  scattered  field  is  determined  by  the  boundary  conditions  on 


the  perfectly  conducting  sphere: 


Btot  -  BinC  +  BSC  = 
r  r  r 

Ekot  -  e£"c  .  E*°  - 


E*nC  +  e!C 


0  ,  r  *  b  <  a 


(5.1) 


13. 


The  scattered  electric  and  magnetic  fields  can  be  derived  from  a  pair  of 
Debye  potentials  u  and  v  which  have  outgoing-wave  eigenfunction 
expansions  similar  to  (k.k)  with  unknown  coefficients.  These  coefficients 
are  determined  by  putting  u*nc,  v*nc  and  vsc  into  (5.1).  In  this  way  we 
easily  get 


sc 

u 


h[l)(kma) 


W> 


hi1)(k»r) 


P^(cos 


e) 


im0-  iu)  t 
m 


sc 


_  00  4  i  i\.  .  ik 

'  Ji  j-t  T7J*I)  ^  (0) 


X  [kah^tka)]' 
mi  m 


[k  bj  (k  b)]’  ,  V 

- - TfT -  hi  (kmr)  ^cose) 

[kb  h1J(kb)]'  11  B  l 

m  i  m 


im0  -  iu  t 


x  e 


m 


,  r  >  b  (5.2) 


These  solutions  are  exact. 


6.  The  Watson  Transformation 

If  we  have  in  mind  the  application  of  the  solution  of  this  problem 
to  the  case  of  a  dipole  antenna  travelling  around  a  planet,  the  expres¬ 
sions  for  the  Debye  potentials  given  by  (5*2)  are  practically  useless. 
The  reason  is  that  in  this  situation  kQb  >>  1  ,  and  the  convergence  of 
the  series  is  extremely  slow.  To  get  a  good  estimate  of  the  sums  we  have 


to  include  terms  up  to  values  of  l  of  the  order  of  kb  .  An  alterna- 

o 

tive  to  direct  summation  is  the  method  of  Watson  transformation^. 

Let  us  first  consider  the  total  Debye  potential  of  the  magnetic 
wave  for  r  >  a  : 

o 

tot  inc  ,  sc 
u  *  u  +  u 


(k  a) 
m 


hi1)(V> 


hi1)(V)  pt(co3 


0)e 


im0-  iw  t 
in 


(6.1) 


Hie  order  of  summation  has  been  inverted  so  that  m  is  now  summed  from 
-•»  to  «  .  This  is  permissible  since  P^(cos  0)  *  0,  ( £—  |m  | )  t  P^(0)  is 
finite  for  |m|  >  l  .  Using  the  addition  theorem 

P  (cos  0  cos  9’  +  sin  6  sin  O'  cos  0) 

Jv 


l 

l 

xm-l 


ifcWl 

U+  m  ) 


p“(cos 


0) 


P^(cos 


0') 


im0 


(6.2) 


we  get 

*«•  •» 

ut0t  can  now  be  rewritten  as 


2tt 

—  |  PJl(sin  0  cos  0,)e‘L“^  d0' 

0 


(6.3) 


tot 

u 


b  l  "lkm 


im(0  -  0' )  -  iw  t 
m 
e 


(6.4) 


where 


15. 


m 


r  21+1  L(2) 

ik  ‘^l  V 


H<2>  (kb) 
t+1  m  (2.) 

(kma)  -  „'(i)  r  ;7  <4  (v> 

(1) 


ht  <kBr)  pt(sin  8  cos  0,)(6.5) 


We  define 


cos  =  sin  Q  cos  0' 


(6.6) 


v  =  £  +  | 


(6.7) 


Then,  since  P^(z)  *  (-1)  P^-z)  for  *  6:1  integer,  Sm  becomes 


m 


2  l  -^T[H(2)(k  a)  -  \T(kmb)  H(l) (k  a)]  h(l)  (k  r) 
v*3/2  V2-  r-  L  V  m  H  1  (k  b)  v  m  J  v-i  m 


Hi2)(kmb) 

■FV1 

v-i 

(-1)  P  1  ( — cos  |J») 
V‘2 


(6.8) 


If  each  term  is  considered  as  a  function  of  the  complex  variable  v  this 
sum  can  be  converted  into  a  contour  integral  in  the  complex  v-plane  (see 
Fig.  1): 

w(l>. 


H(2)(k  b) 


”  l  v2-  J  l  V  B  H.(1,(k.b)  v  ”  J 


(k  r) 


v  m 


1 

i-  £  m 


cos  vit 


P  .  (-cos  \J>)  dv 
v-i 
2 


(6.9) 


The  contour  C  encloses  the  positive  real  axis  along  which  are  the  posi¬ 
tive  poles  of  1/cos  vit  .  All  other  singularities  of  the  integrand  are 
to  be  excluded.  The  residues  at  these  poles  give  back  the  sum  (6.8), 


We  see  that  the  factor  involving  the  Hankel  functions  is  also  an  even 
function  of  v  .  Thus  the  entire  integrand  is  an  odd  function  of  v  .  In 
this  case  the  integration  along  the  lover  branch  of  the  contour  C  can 
be  replaced  by  one  along  its  image  vith  respect  to  the  origin  in  the  second 
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quadrant,  as  shown  by  the  broken  line  in  Fig.  1.  The  integral  is  now 
to  be  evaluated  along  a  contour  lying  Just  above  the  real  axis: 


H  Vd) 

s*  ’  -l  1  -tl  [Kv2)  V>  -  -rry 


ci  ’  f 


H(2)(k  b)  (  v 


u  (kmb) 
v  m 


h(l),(k  r) 
m 


COS  VTT  V-  — 


P  1  (-cos  )dv 


(6.1U) 


In  a  similar  way  we  consider  the  total  Debye  potential  of  the  electric 
wave  for  r  >  a: 


tot 


inc  ,  sc 
v  +  v 


&  i  _  ilvlmli!  ^'(0)~r 


"o  m=-®  t»l 


[km* 


tkJ>  J.<kJ>M  (!) 


m 

[k  b  hf’fk  b)]' 
m  *.  m 


(v  hr  <kma»'; 


* 


im0-iu  t 

x  P^(cos  Q)  e  (6.15) 


Using  the  addition  theorem  (6.2)  we  get 


2tt 


cv-  \  snsbr  *’ 

0 


(6.16) 


Thus 


2ir 


tot 


P_ 

♦ire 


\  *'  fe  l  -f 

L  m«-« 


ik  im(0-0’ )-  iu  t 
m  m 


8  3(cos  iji)  im 


T  cos  0 


(6.17) 
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where 


T  -  2  [  f[^T7  H(2)(k  a)]'  - 

m  v-3/2  v2-  ±  L  m  m 


[  /FT  H^fk  b)]' 

m  v  m 

[  /FT  H^fk  b)]1 

m  v  m 


[  /FT  H^(k  a)] 
m  v  m 


■]  ■:,a  >*■ 


r)(-l)  (-cos  if/) 

v-  g 

(6.18) 


T  can  be  converted  into  a  contour  integral  by  the  Watson  trans format ion: 
in 


t  v  r  _  (c>)  [/ft H'2)(k b)]' 

P  ■  -i  ”5 — r  [ /k  a  H  2^(k  a)]' - - - /-.■') -  [/k"i 

m  Jv2_  11  m  v  m  [/FT  Hp:)(kinb)]' 


a  H(l)(k  a]' 
re. 


h(1)  (k  r) 
v-i  " 


COS  V1T  V  -  2 


P  ,  (-cos  ij/)  dv 


(6.19) 


The  contribution  from  the  double  pole  at  v  *  g  be  subtracted  off  is 


.  f  r  ,0.  EiTbH^^kb)]'  ... 

r  *  2  rrl — 2-r  [/Ta  H  2  (k  a)]' - -----  -frx — 2 -  [/Ta  O  (kma)]' 

m  3v|v+H  m  v  »  [vFT  H(1)(k  b)]'  m  v  m  J 

v  c.  m  v  m 


h  (k  r)  p  _  (—cos  if/)) 

m  v-a 


(6.20) 


The  integrand  in  (6.19)  is  again  an  odd  function  of  v  .  The  integral 
can  be  evaluated  along  a  contour  lying  just  above  the  reed  axis: 
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m 


“1 


v2- 


l 

IT  L 


H^(k  a)J 
in  v  m 


[✓Tb  H^2)(k  b)]' 

a  v  m 

[ZkT  H^tkb)]1 
m  v  m 


X 


X 


(k  a)]’ 
a 


h(l),(ltmr) 
V~? 
cos  vir 


(-cos  ip)  dv 


(6.21) 


We  now  consider  the  contribution  from  the  pole  at  v  ■  1/2  .  After 
some  straightforward  calculations  we  find  that  the  parts  of  utot  and 
vtot  due  to  om  and  Tffi  Jointly  give,  upon  substitution  in  ( U . 3 ) »  iden¬ 
tically  zero  contribution  to  the  total  field.  From  now  on  these  parts 
will  be  discarded. 

7.  Separation  of  the  Geometrical  Optics  Field 

We  want  to  evaluate  the  contour  integrals  (6.14)  and  (6.21)  asymp¬ 
totically  in  the  limit  of  a  very  large  sphere:  kQb  >  >  1  .  For  this 

2 

purpose  it  is  most  elegant  to  use  a  method  due  to  Franz  to  separate  the 
total  field  into  two  parts.  One  part  can  be  identified  with  the  geomet¬ 
rical  optics  field,  and  the  other  is  in  the  form  of  a  series  of  damped 
"creeping  waves". 

As  will  be  seen  below,  for  kb  >>  1  ,  most  of  the  contributions 

m 

to  the  integrals  (6. lb)  and  (6.21)  come  from  large  values  of  v  of  the 
order  of  kb.  From  the  identity 


Pv(— cos  ip)  ■  eivir  Pv(cos  ip)  -i  sin  vw[Pv(cos  ip)-i  ~  Qv(cos  ^)]  (7.1) 


and  the  asymptotic  formulas 
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P  (cos  ip)  -V  /  2  "  cos [  ( v  +  -  r-] 

v  V  irv  sin  t|/  2  k 


f  Q41(cos  $)  *  J —  -  cos [  ( v  +  +  r*]  ,  |v  sin  (|>|  >>  1  (7.2) 

V  irv  sin  0  4 


IT  V 


we  get 


P  (-cos  1 1>)  +  C— ■—  e  *  i^+e1(V  2>P  (cos  i|») 
v  -  *■  /2ir  v  sin  i|)  v-^- 


|v  sin  <J/|  >>  1  (7-3) 

Thus  for  sin  i|»  ^  0  we  can  substitute  (7-3)  into  the  integrals,  each  of 
which  is  now  separated  into  two  parts: 


m 


S8  ♦  scr 

m  a 


'  T 


m 


T8  +  TCr 
m  m 


(7.U) 


where 


3»  ■/ 


_  j  n 

2 - ilT 

e 


iTsin  ♦ 


,(«(v).&£i  .(Du 

.  v-r'-'1  m  nu)  v  m  . 

4  v  m 


f  JL-  ’ 

J  2  1  J 

i  V  -  T  L 


4 

m  y  i 


x  h^  (k  r)  eiv<^dv 
1  m 

V"2 

C/kT  H(2)(kb)l' 


it  sin)p 


^  V2-  [k_b  H  1  (k  b)]' 


C1  '  * 


[£TH(l)(ka)]'  h(l\  (k: 
m  m  i  m 

v-  *■ 


m  v  m 


r)  elvl"  dv 


(7.5) 
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and 


Scr 

m 


“I  v2-JH‘ 


pN  H[2 !(ltb)  ,  » 

2)(k  a)  -  irr-H^ika) 

m  H  1  (k  b)  v  m 
v  m 


h(l)  (k  r) 
v-±  m 

x  - 2 -  p  ^ ( cos  tp)  e*V7T  dv 

cos  vtt  v-  — 


rcr 

m 


i  v  -  rl 


[^Tb  H(2)(k  b ) ] ' 
m  v  m 

[JTb "Hp^{k  b)]' 

m  v  m 


m  v  m 


h(l)l(k mr) 


i  VTT 


_  P  ^cos  )  e  dv  .  (7*6) 

cos  vir  v”  2 


In  the  next  section  it  will  be  seen  that  Sg  and  Tg  in  (7.5)  yield  the 

m  m 

geometrical  optics  field. 

It  be  shown  that  the  integrands  in  (7-6)  tend  to  zero  at  infinity 

in  the  upper  v-plane  for  r  >  a  >  b  ,  except  in  the  neighborhood  of  the 

zeros  of  H^^(k  b)  or  [/k  b  H^(k  b)]'.^  For  k  b  >>  1  these  are 
v  m  m  v  m  m 

located  approximately  at 

vs  *  V  +  Vkmb)1/3  ei"/3  •  (7.7) 

where  the  A’s  are  real  positive  constants  which  are  different  for  the 
magnetic  and  electric  waves.  We  observe  that  these  zeros  lie  in  a  row 
in  the  first  quadrant  (see  Fig.  2).  They  are  simple  poles  of  the  inte¬ 
grands  and  are  the  only  singularities  in  the  upper  v-plane.  The  straight 
contour  can  now  be  deformed  to  the  contour  C2  which  encloses  these 
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poles,  as  shown  in  Fig.  2.  The  residues  obtained  from  these  poles  form 
an  infinite  series  of  "residue  waves"  which  are  rapidly  convergent  for 
large  k^b  .  Each  term  of  the  series  can  be  interpreted  as  an  exponen¬ 
tially  damped  "creeping  wave".  It  must  be  remarked  that  we  have  tacitly 
assumed  that  the  main  contributions  to  utot  and  vtot  come  from  terms 

with  small  m's  such  that  k  b  =  k  b  +  mJCb  -  k  b  >>  1  . 

mo  o 

In  what  follows  we  will  confine  our  attention  to  the  geometrical 

optics  field,  since  it  is  more  interesting,  besides  being  the  dominant  part 

for  k  b  >>  1  . 
o 

8.  The  Geometrical  Optics  Field 

We  first  consider  the  integral  Sg  in  (7-5).  We  will  evaluate  it 

m 

by  the  method  of  stationary  phase.  This  integral  in  general  has  one  or 
two  stationary  points  on  the  real  axis  in  the  neighborhood  of  v  *  kmb  . 

In  order  to  simplify  the  problem  we  take  the  limit  r  >>  a  so  that  we  will 

er 

be  calculating  the  far  field.  It  is  now  permissible  to  apply  to  S°  the 
Hankel  asymptotic  formula 

(l)  ,  i(kr-vf-f) 

h  1  i(k  r)  *  —  e  2  U  kmr  »  v  (8.1) 

v  -  *■  m 


The  case  where  r  is  of  the  same  order  as  a  should  present  no  special 
difficulty;  but  the  calculations  that  follow  would  be  more  cumbersome.  Sub¬ 
stituting  (8.1)  into  Sg  we  get 


3g  «/ 
m  v  i 


m 


ik  r  -  i^-  • 


w  sin  ij; 


» _ i  rjGL 

k  r  I  2  1 
m  1  v  -  r 


„<*>»  a) 

v  m  H  1  (k  b)  v  m 
v  m 

iv(*  -  f) 
x  e  dv 


(8.2) 
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To  locate  the  stationary  points  we  apply  to  (8.2)  the  Debye  asymp¬ 
totic  formulas  for  the  Hankel  functions  in  the  limit  when  v  and  z 
are  large  and  comparable: 


H(l)(z) 

v 

h(2)(z> 


)'  Art 


'V  +i 


,/£ _ i- 

/  ’  (v2-  z2: 


/  2  2  -1  v  TT. 

+  i(vz  -v  -  v  cos  —  -  f-) 

e  2  ,  v  <  z 


v+  /?-  z2  /  2  2 


-  ✓  V  -  z 


,  V  >  z 
(8.3) 


with  the  condition 


„  ~1  V  TT 

0  <  cos  —  <  x 
Z  £ 


(8.U) 


We  split  SK  into  three  integrals 
m 


ik  r  -  i£ 

- - - -  m  p 

s?  =,/■  ■  5 — r-r —  [In  +  Io  +  Ij 


m  u  it  sin  tlh  k  r 

*  m 


1  2  3 


(8.5) 


where 


i  v  -  r 


iv(*  -  J-) 


(8.6) 


7  *^7  \  (kmb>  (1)  iv(*  -  §■) 

-  ~pi  fl)  E  e  < 

kJb  v  -  ¥  Hi  (kmb) 


(8.7) 


H(2)(k  b) 


_  f  ^  g(i)(lt  a) 

J  v2-f7^7  v  (  m  ) 


ivU  -  f) 

e  dv 


v  m 


(8.8) 


Substituting  the  appropriate  Debye  asymptotic  formulas  in  these  integrals 


we  get 


2k, 


!  mfl  ii  7  AT  1  1 

1  »*  _J  v2-  $  0c=a2-  v2)1A  ' 


l-{i 


,  2  2  ?  .  -1  v  .  /  ,  it  v , 

kma  -v  ♦  v  cos  —  +v(*-^)] 


dv 


(8.9) 


/_  09 

2  *  f  1 

ir*  I  2  1  7^2  1717^ 

k  b  v  -  IT  (kma  -  v  > 

m 


i[»42a2-  v2-  v  cos-1  j^+  v(<h  |-)  ] 
e  m  dv 


(8.10) 


k  b 
.  JL  m 

ik 


fi  'ii  f  _ l 

V"  '  J  .2  I  ,.2.2  .2,1/11 


c  j.  ,,  c.  c.  cij 

v  -  (k_a  -  v  ) 


m 


x  e 


,  r  f.2  2  2  n.f.  2.  2  2  -I  v  _  —1  V  [  |  Hii 

ill/k  a  -  v  -  2yk  b  -  v  -  v  cos  - - +  2v  cos  r-rr  +  v(ij>  -  tt)J 

m  '  m  ka  kb  2 

m  m 

i  dv 


k  a  >  k  b  >>  1  .  (8.11) 

m  m 


In  writing  down  the  above  expressions  we  have  used  the  fact  that  1^  and 
1^  have  no  stationary  points  to  the  right  of  k^a  . 

Each  of  the  three  integrals  has  at  most  one  stationary  point.  For 
the  integral  I1  the  location  of  its  stationary  point  is  given  by  the  solu¬ 
tion  of  the  trigonometric  equation 

005-1  iTa  +  ♦  “  2  =  0  (8.12) 

mB 

Because  of  the  condition  (8.4),  this  equation  has  a  solution  only  when 
<  j  .  When  this  is  the  case  we  have 

v  ■  k  a  sin 


(8.13) 


Af.t.mmwnr  *  .  9«wr 
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Thu9  evaluating  1^  by  the  well-known  method  of  stationary  phase,  we 


get 


I,  « 


-ik  a  cos  ip 
m 


(k  a  sin  ip) 
in 


3/2 


♦<5 


<P  > 


(8.1U) 


In  the  same  way  the  location  of  the  stationary  point  of  1 2  is  given  by 

-1  v 


-  cos  ~  kT  +  *  -  I 

m 


(8.15) 


Because  of  (8.M,  this  equation  has  a  solution  only  when  ip  >  —  .  Then 
the  stationary  point  is  given  by 


v  *  k  a  sin  ip 
m 


(8.16) 


This  must  be  greater  than  kb  or  else  the  point  lies  outside  the  lower 

m 

limit  of  integration.  Hence  we  have  the  additional  condition  for  the 
existence  of  the  stationary  point: 


sin  ip  >  -- 

a 


(8.17) 


Therefore  we  get 


(k  a  sin  ip) 
m 


»  <  2 


-ik  .  cos  * 

57s  e  .  a  <  *  <  *  -  .1-  X 


,  it  -  sin-1  —  <  ip  <  it 
*  a 


(8.18) 
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The  sum  1^  +  1^  has  a  very  obvious  interpretation  in  geometrical 

optics.  It  represents  the  incident  wave.  ^  as  defined  in  (6.6)  is  the 

angle  between  the  directions  (0,0')  and  (~,0)  .  Suppose  a  point  source 

of  frequency  is  placed  at  the  point  (a,j,Q)  as  shown  in  Fig.  3.  The 

geometrical  shadow  region  is  a  cone  with  apex  at  the  source  and  angle 

sin-1  —  .  Inside  the  shadow  region  I,  +  I„  is  zero.  Outside  this 
a  i  d 

region  I,+  I„  contributes  a  term  to  S®  proportional  to 
j.  d  m 

,  ik  (r  -  a  cos  iii) 

1  m  r 


which  indeed  describes  the  far  field  of  a  point  source  displaced  from  the 
origin  by  a  distance  a  . 

We  now  want  to  show  that  I^  gives  rise  to  the  reflected  wave. 

This  integral  has  a  stationary  point  given  by 


-1  v 


-1  v 


*  coe  .  kT  ♦  2  003  iTb  +  *  -  2 

m  m 


(8.19) 


We  denote  the  solution  by 


v  *  kb  sin  y  ,  0  <  y  <  £ 

m  c 


(8.20) 


Then  Y  is  the  angle  of  reflection  for  a  ray  scattered  into  a  direction 
inclined  at  an  angle  ip  with  the  x-axis.  This  can  readily  be  seen  from 
the  ray  diagram,  Fig.  3.  Without  loss  of  generality  we  have  taken  the 
z-x  plane  to  be  the  plane  of  reflection.  From  this  diagram  we  get 


b  sin  y  *  a  sin  a 


(8.21) 


where  a  is  the  angle  the  incident  ray  makes  with  the  x-axis.  Using  this 
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result  in  (8.19)  we  get 

a  -  2y  +  ii  =  0 


(8.22) 


which  is  just  the  law  of  reflection.  Eliminating  a  from  (8.21)  and 
(8.22)  we  obtain 


sin(2y  -  iM  _  b 
sin  y  a 


(8.23) 


This  is  a  quartic  equation  in  sin  Y  which  can  always  be  solved.  The 
exact  expression  for  the  solution,  however,  is  too  complicated.  For 
Y  >  ^  ,  that  is,  inside  the  geometrical  shadow  region  (8.19)  has  no  solu¬ 
tion  and  there  is  no  stationary  point.  Thus  we  have 


b  cos  y 


(k-b  sin  '  '1 2  V^-"b2sin2Y  -  b  cos  Y 


m 


x  e 


ik^/a2-  b2sin2Y  -  2b  cos  y) 


_  ,  -1  b 

,  0  <  ip  <  tt  -  sin  — 
Y  a 


0 


.-lb 

tt  -  sin  —  <  0)  <  ti 
a  r 


(8.2U) 


Summarizing  the  foregoing  results  we  have 
ik  r  _  r 

sg  =  -i  sJL  r I  1 

m  r  V  n  sin  \p  5/2 


-ik  a  cos  1 1 ) 
m 


Yr'  L (a  sin  \J») 
m  r 


3/2 


(b  sin 


\>3'2  i: 


b  cos  Y 


ik  A&  b2sin2y  -  2b  cos  Y ) 
m 


2ya  -  b  sin  y  -  b  cos  Y 


0  <  <  tt  -  sin’^  — 

a 


,j mmn 


o 


v  -  sin”^  —  <  <J<  <  ti 
a 


(8.25) 


For  the  electric  wave  we  get  in  a  similar  manner 


IK  r  t 

m  / - q - i/k  a  f  -i  cos  iji  -ik  a  cos  ^ 

mS  u  i  e  /  2  —  m  - . —  -  111  + 

“  1  r  /.  sin*  fc5/2  [(a  ,ln  ^ 

m 

+  t  vggHS  /  b  an.  v  ^ikn(/a2-bSsi„2  -2b  cosy) 

a(b  sin  y)3^  y  g |/a2-b2sin2Y  -  b  cos  y 

0  <  i/j  <  ti  -  sin"'*'  — 

a 


.  -1  b 

ti  -  sin  —  <  il»  <  ir 
a 


(8.26) 


From  this  we  get 


3 (cos  ^  “  -i 


ik  r 

e  m  rr 


a  -k  a  cos  ^  -ik  a  cos  il> 
T  m  a  m 

l^Lu  sin  *>3/2 


r-  y  s  sin  *  k5/2L(a  sin  *>3/2 

m 


k  42-  b2sin2Y 
m _ 


h  cos 


ik  b2sin2  -  2b  cos  Y  )j 


a  sini|>  (b  sin  y)1^2  v  2  /a2-  b2sin2Y  -  b  cos  Y 


o  .  ,  -lb 

0  <  ^  <  ir  -  sin  —  , 

a 


0  , 


n  -  sin”^  —  <  <  ti  . 

a 


(8.27) 


In  the  above  we  retained  only  the  highest  order  terms  in  k^a  . 
Substituting  (8.25)  into  (6.U)  we  get  for  the  incident  part 


4 tic  ra* 


f  d£  1  r  m  1LV"VV 

~~2  p7  L  iT  6 
■*  sin  w  m*-«  m 


i[k  r-u  t-k  a  cos  i|»  +  m(0  -  0' )] 


f  *  IJt' 
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where  the  integral  is  taken  only  over  the  illuminated  region,  that  is, 
over  the  range 


-cos 


.  J~ 2  v2 

-1  -fa  -  b 
a  sin  9 


<  0'  <  cos 


-1  -  l/a^-  b2 


a  sin  9 


From  (8.28)  we  get 


3  -.inc  _  P _ 3  I  Jrtl,  e 


i(k  r  -  w  t  -  k  a  cos  i|>) 
o  o  o 


T,77^ 

o 


d0* 


sin2<J> 


.  «  im(Kr  -  Qt  -  Ka  cos  \p  +  0  -  0'  ) 

*  T7  I  • 


IJlS-00 


(8.29) 


By  (2. U )  we  have 


1 

2  r 


•  im(Kr  -  Qt  -  Ka  cos  ^  +  0  -  0'  ) 


I  • 

m=-<>“ 


6(Kr  -  Qt  -  Ka  cos  ip  +  0  -  0' ) 


(8.30) 


The  position  of  the  peak  of  the  delta- function  i3  given  by  the  solution  of 
the  equation 

Kr  -  Qt  -  Ka  sin  9  cos  0'  +  0  -  0*  =  0  (8.31) 

In  most  situations  Ka  =  ~  <<  1  .  (8.31)  can  be  solved  by  iteration: 


0 '  =  0  +  Kr  -  Qt  -  Ka  sin  9  cos(0  +  Kr  -  Qt)  + 


(8.32) 


In  what  follows,  for  simplicity,  we  will  ignore  all  terms  in  Ka.  Therefore, 
(8.29)  becomes 


30. 


i[k  r-(o  t-k  a  sin  0  cos(0  +  Kr-  fit)  ] 
3  inc  P  3  e  0  0 

■  ■  Y\X  B  .  ■■■  —  — i  i  ■■■  ■  i  ■  i  i  ■  mmm  — ■—  ■■■ 

Une  a2  ^  1  -  sin2©  cos2(0  +  Kr  -  £>t ) 

o  *- 


(8.33) 


when  0  ♦  Kr  -  fit  lies  in  the  illuminated  region  and  zero  otherwise. 
(8.33)  is  equivalent  to 


P _ sin  0  sin(0  +  Kr-fit) 

l*ne  a  1  -  sin2©  cos2(0  +  Kr  -  fit ) 
o 


x  e 


i[k  r-u>  t-k  a  sin  ©  cos(0  +  Kr-  fit )  ] 
o  o  o 


(8.3*0 


Carrying  out  the  same  calculations  for  the  reflected  part  we  get 


rure  ■  -  P  sin  ©  sin(0  +  Kr  -  At)  / a  sin  b  cos  y 

^’eoa  1-  sin2©  cos2(0 +  Kr- fit)  Vb  sin  Y  2\/a2-  b2sin?Y  -  b  cos  Y 


x  e 


i[k  r-wt+k  (\jaf-  b  sin  y-2b  cos  y)]  (8.35) 

o  o  o 


Hie  angles  Y  and  are  to  be  evaluated  at  0'  =  0  +  Kr  -  fit  .  The 
whole  expression  is  zero  in  the  shadow  region.  Similarly  for  the  electric 
wave  we  get 


P  cos  0  sin  0  cos(0  +  Kr  -  fit) 

^£oa  1-  sin2©  cos2 (0  +  Kr  -  fit) 


i[kr-wt-ka  sin  ©  cos  (0  +  Kr  -  fit )  ]  (8.36) 

ooo 
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ref  P 


/  b2  2 
cos  0  |/1-  —  sin  Y 

_ *' _ a  d 


a  sin 


b  cos  Y 


+7T€oa  1-sin2©  cos2(0+K r-  ftt)  V  b  sin  Y  2{I2-b28ln2Y  - 


b  cot  Y 


i[k  r-o)  t+k  ( ]/a  -  b  sin^Y  -  2b  cos  y)] 

«  «  0  0  1  (8.37) 


>.  Nature  of  the  Solution 


We  recall  that  the  solution  is  zero  when  the  direction  (0,0  +  Kr-flt) 
lies  inside  the  shadow  cone  of  Fig.  3.  0  +  Kr  -  fit  is  Just  the  difference 
between  the  azimuthal  angle  0  of  the  direction  of  observation  (0,0)  and 
the  retarded  azimuthal  angle  H(t-“)  of  the  revolving  dipole.  Thus  we 
have*  the  simple  conclusion  that  the  total  field  is  zero  when  the  observa¬ 
tion  point  lies  inside  the  retarded  shadow  region  of  the  dipole. 

We  substitute  ruinc  and  rvinc  in  (8.31*)  and  (8.36)  into  (!».i)  and 
obtain  the  asymptotic  form  of  the  incident  field  in  the  illuminated  region: 


c  i[kr-u>t-ka  sin0  cos(0+  Kr  -  ftt)] 

O  i  a  0  0  0 

——  sin  0  e 


4ire  r 
o 


(9.1) 


When  fl  ■  0  these  field  components  coinc-'de  with  those  of  a  Hertz  dipole 
fixed  at  the  point  x*a,  y*0,z=0  and  oriented  parallel  to  the 
z-axis.  We  can  define  an  instantaneous  frequency  of  the  incident  wave  by 
differentiating  the  exponential  with  respect  to  t  : 


to  ■  u  [l  +  J  sin  0  sin(ftt  -  Kr  -  0)] 
o 


(9.2) 


where 


32. 


i 


(9.3) 


Here  fia  is  Just  the  velocity  of  the  dipole,  and  sin  0  sin^fit  -  Kr  -  0) 

is  the  cosine  of  the  angle  between  the  velocity  of  the  dipole  and  the 

direction  of  observation  at  the  retarded  time  t  -  —  .  (9*2)  therefore 

c 

agrees  with  the  Doppler  formula  for  the  frequency  shift  of  a  moving  source. 

Similarly  we  calculate  the  far  field  of  the  reflected  wave  in  the 
illuminated  region  from  ruref  and  rvref  in  (8.35)  and  (8.37): 


0 


,re 


2  _ 

f  „ref  P  o  /b  sin  y  b 

’cB«  ^7 

sin  9  sin20’  +  cos29  sin  0'Jl  -  ^  sin£ 


21 

sin  Y  -  b  cos  y 


x  e 


sin2'!' 


i[k  r-w  t  +  k  (  a2-  b2sin2Y  -  2b  cos  y)] 
o  o  o 


,ref 


„ref  pko  /b  sin  Y  b  cos  Y 

-°B0  ■: -  1/  a  sin  ->  "" "  * - - 

iiir  p  r  * 


Uirc  r 
0 
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b  sin  Y  -  b  cos  Y 


cos  9  sin  0'(sin  9  cos  0' 


-P? 

v  a 


sin2Y ) 


sin2^ 


i[k  r- u)  t  +  k  (/a2-  b2sin2Y  -  2b  cos  Y)]  (9«*0 

o  o  o 
e 
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33. 


The  angles  <p,  y  and  0'  are  defined  in  terms  of  r,  0,  0  and  t  by 
(6.6),  (8.23)  and  (8.32)  respectively.  The  instantaneous  frequency  of 
the  reflected  wave  is  given  by 

</ef  =  oi  [l  +  B'sin  0  sin(ftt  -  Kr  -  0)] 
o 

where 

>  fib  sin  Y 
B  c  sin  ip 

Thus  the  reflected  wave  is  not  of  the  same  frequency  as  the  incident  wave. 
Unlike  B  in  ( 9 • 3 )  6'  is  a  periodic  function  of  time.  The  quantity 
b  sin  y/sin  is  the  intercept  of  the  ray  -eflected  into  the  direction 
(0,0)  when  projected  backwards,  on  the  retarded  radius  vector  of  the  moving 
dipole  as  can  be  seen  from  Fig.  3.  Despite  the  apparent  similarity  of  (9.2) 
and  (9*5)  it  does  not  seem  to  be  possible  to  interpret  the  reflected  wave 
as  due  to  a  virtual  moving  point  source. 
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